Determining optimal gate sizes by using a numerical solver

ABSTRACT

Systems and techniques are described for optimizing a circuit design by using a numerical solver. The gates sizes are optimized by modeling a set of gate optimization problems and solving the set of gate optimization problems by using a numerical solver. Modeling each gate optimization problem can include selecting a portion of the circuit design that includes a driver gate that drives one or more inputs of each gate in a set of gates, and modeling a gate optimization problem for the portion of the circuit design based on circuit information for the portion of the circuit design. A differentiable objective function for delay can be created using numerical models for the delays in the circuit. In some embodiments, gradients of the differentiable objective function can be computed to enable the use of a conjugate-gradient-based numerical solver.

TECHNICAL FIELD

This disclosure relates to electronic design automation (EDA).

More specifically, this disclosure relates to determining optimal gate sizes by using a numerical solver.

BACKGROUND Related Art

The goal of circuit synthesis is to convert a high-level description of a circuit design into an implementation that meets a set of timing constraints, and at the same time optionally optimizes one or more metrics, such as area, leakage power, etc.

Some circuit synthesis approaches create an initial circuit design at a given abstraction level (e.g., a logical or physical design). Next, a cell is identified in the circuit design for optimization based on the metrics that are desired to be optimized. An optimal size for the identified cell is then determined by iteratively replacing the identified cell with functionally equivalent cells that have different sizes (this optimization process is also referred to as “sizing the cell,” “sizing the gate,” etc.). For each replacement cell size that is tried, the circuit synthesis approach updates timing information, and rejects cell sizes for which one or more timing constraints are violated. The iterative optimization process typically terminates after the optimization process has executed for a certain number of iterations or for a certain amount of time. One aspect that further makes this process computationally expensive is that modern technology libraries have many sizes associated with each gate type with carrying driving capability and/or multiple threshold voltages.

Unfortunately, such iterative trial-and-error based circuit synthesis approaches either take too long to complete and/or produce poor quality results for large circuit designs in which timing constraints are checked across many process corners and modes.

The cell delay can be represented using the following expression: d=R·C _(o) +p,  (1) where, R is the output resistance of the cell, C_(o) is the output loading, and p is the parasitic delay of the cell. Equation (1) can then be rewritten as:

$\begin{matrix} {{d = {{\left( {R \cdot C_{i}} \right) \cdot \left( \frac{C_{o}}{C_{i}} \right)} + p}},} & (2) \end{matrix}$ where, C_(i) is the input capacitance presented by the cell at one of its input terminals.

A circuit synthesis approach that is different from iterative trial-and-error approaches is described in Ivan Sutherland, Robert F. Sproull, and David Harris, Logical Effort: Designing Fast CMOS Circuits, Morgan Kaufmann, 1999. This approach uses the following linear cell delay model: d=g·h+p,  (3) where, g represents the logical effort, h represents the electrical effort, and p represents the parasitic delay of the cell. The logical effort captures the effect of the cell's topology on its ability to produce output current. The logical effort is independent of the size of the transistors in the circuit. The electrical effort describes how the electrical environment of the cell affects performance, and how the size of the transistors in the cell determines its load-driving capability. The parasitic delay is a form of delay overhead that accompanies any gate. The logical effort g and the parasitic delay p can be normalized, e.g., they can be normalized with respect to the logical effort and parasitic delay values, respectively, of a reference inverter.

Comparing Equations (2) and (3) we note that the first term (R·C_(i)) represents the logical effort g, and the second term

$\left( \frac{C_{o}}{C_{i}} \right)$ represents the electrical effort h.

The delay model that is used by a circuit synthesis approach must be accurate, i.e., it must accurately model the actual cell delays. If the delay model is inaccurate, the circuit implementation produced by the circuit synthesis approach will likely contain many timing violations. In the above-mentioned work by Sutherland et al., the authors state that the linear delay model is approximate, i.e., it does not accurately model the actual cell delays.

SUMMARY

Some embodiments described herein provide techniques and systems for determining optimal gate sizes by using a numerical solver. The gates sizes are optimized by modeling a set of gate optimization problems and solving the set of gate optimization problems by using a numerical solver (e.g., a conjugate-gradient based numerical solver).

Modeling each gate optimization problem can include selecting a portion of the circuit design that includes a driver gate that drives one or more inputs of each gate in a set of gates, and modeling a gate optimization problem for the portion of the circuit design based on circuit information that can include generic logical effort values of each gate in the set of gates, an input capacitance value and a specific logical effort value of the driver gate, and wire resistance and capacitance values of a net that electrically connects an output of the driver gate with one or more inputs of each gate in the set of gates. The wire resistance value and the wire capacitance value of the net can generally be determined based on any wire load model. Examples of wire load models include, but are not limited to, an Elmore wire load model, a best case wire load model, a worst case wire load model, and a balanced-tree wire load model.

In some embodiments, optimal gate sizes can be estimated based on a set of closed-form expressions, and the estimated gate sizes can be provided to the numerical solver as a seed to solve the gate optimization problem.

In some embodiments, multiple optimization iterations can be performed, wherein optimal gate sizes determined in a given optimization iteration are used to model optimization problems in the next optimization iteration.

BRIEF DESCRIPTION OF THE FIGURES

FIG. 1 illustrates a numerical synthesis process in accordance with some embodiments described herein.

FIG. 2 illustrates a reverse-levelized processing order in accordance with some embodiments described herein.

FIG. 3A how a gate optimization problem can be modeled in accordance with some embodiments described herein.

FIGS. 3B-3D illustrate a best case wire load model, a worst case wire load model, and a balanced-tree wire load model, respectively, in accordance with some embodiments described herein.

FIG. 4 illustrates a process for optimizing a circuit design in accordance with some embodiments described herein.

FIG. 5 illustrates a computer system in accordance with an embodiment of the present invention.

DETAILED DESCRIPTION

The following description is presented to enable any person skilled in the art to make and use the invention, and is provided in the context of a particular application and its requirements. Various modifications to the disclosed embodiments will be readily apparent to those skilled in the art, and the general principles defined herein may be applied to other embodiments and applications without departing from the spirit and scope of the present invention. Thus, the present invention is not limited to the embodiments shown, but is to be accorded the widest scope consistent with the principles and features disclosed herein. In this disclosure, when the term “and/or” is used with a list of entities, it refers to all possible combinations of the list of entities. For example, the phrase “X, Y, and/or Z” covers the following cases: (1) only X; (2) only Y; (3) only Z; (4) X and Y; (5) X and Z; (6) Y and Z; and (7) X, Y, and Z.

Overview of an Electronic Design Automation (EDA) Flow

An EDA flow can be used to create a circuit design. Once the circuit design is finalized, it can undergo fabrication, packaging, and assembly to produce integrated circuit chips. An EDA flow can include multiple steps, and each step can involve using one or more EDA software tools. Some EDA steps and software tools are described below. These examples of EDA steps and software tools are illustrative purposes only and are not intended to limit the embodiments to the forms disclosed.

Some EDA software tools enable circuit designers to describe the functionality that they want to implement. These tools also enable circuit designers to perform what-if planning to refine functionality, check costs, etc. During logic design and functional verification, the HDL (hardware description language), e.g., SystemVerilog, code for modules in the system can be written and the design can be checked for functional accuracy, e.g., the design can be checked to ensure that it produces the correct outputs.

During synthesis and design for test, the HDL code can be translated to a netlist using one or more EDA software tools. Further, the netlist can be optimized for the target technology, and tests can be designed and implemented to check the finished chips. During netlist verification, the netlist can be checked for compliance with timing constraints and for correspondence with the HDL code.

During design planning, an overall floorplan for the chip can be constructed and analyzed for timing and top-level routing. During physical implementation, circuit elements can be positioned in the layout (placement) and can be electrically coupled (routing).

During analysis and extraction, the circuit's functionality can be verified at a transistor level and parasitics can be extracted. During physical verification, the design can be checked to ensure correctness for manufacturing, electrical issues, lithographic issues, and circuitry.

During resolution enhancement, geometric manipulations can be performed on the layout to improve manufacturability of the design. During mask data preparation, the design can be “taped-out” to produce masks which are used during fabrication.

Numerical Synthesis

FIG. 1 illustrates a numerical synthesis process in accordance with some embodiments described herein.

Technology library 102 includes library cells and discretized non-linear delay models for the library cells. The library cells in the technology library 102 can be used to create and/or optimize a circuit design. The term “library cell” refers to a cell in a technology library. The term “library cell” is to be distinguished from the term “cell instance” which is an instantiation of a library cell in a circuit design. In some embodiments, a discretized delay model models the delay for a timing arc of a library cell.

The term “library cell type” refers to a particular logical functionality. Examples of functionalities include, but are not limited to, “AND,” “OR,” “XOR,” “multiplexor,” etc. The term “library cell” refers to a particular implementation of a library cell type. A technology library can be viewed as a collection of library cells of one or more library cell types. For example, a technology library may include multiple sizes of an inverter. In this example, the term “inverter” is the library cell type, and the term “library cell” refers to an inverter implementation of a particular size.

A numerical delay model is a delay model that can be used by a numerical solver to optimize a cost function. A linear delay model is one example of a numerical delay model. More generally, an analytical formula that represents the delay behavior is an example of a numerical delay model.

The term “specific numerical delay model” refers to a numerical delay model that models the delay behavior of a particular timing arc of a particular library cell, or the delay behavior of a particular library cell. A specific numerical delay model is defined using a specific logical effort value and a specific parasitic delay value.

The term “generic numerical delay model” refers to a numerical delay model that models an aggregated delay behavior of either a timing arc of a library cell type or the library cell type itself. A generic numerical delay model is defined using a generic logical effort value and a generic parasitic delay value.

Examples of systems and techniques that can be used to determine specific and generic numerical delay models (which can be used to perform operation 104) are described in U.S. patent application Ser. No. 13/450,178, entitled “Numerical delay model for a technology library cell and/or a technology library cell type,” the contents of which are herein incorporated by reference in their entirety.

Some embodiments determine numerical delay models based on technology library 102 (operation 104). Output 106 of operation 104 can include a set of library cells and a specific numerical delay model for each library cell in the set of library cells and a generic numerical delay model for each library cell type. The set of library cells in output 106 typically includes all library cells in technology library 102, but in some embodiments certain library cells may be removed if desired.

Certain library cells can then be excluded (operation 108) to obtain a subset of library cells. Each library cell in the subset of library cells has a specific numerical delay model that was computed in operation 104. In some embodiments, a new generic numerical delay model can be determined for each cell type based on the specific numerical delay models corresponding to the subset of library cells. These generic numerical delay models are “new” because they are based on the subset of library cells as opposed to being based on the set of library cells that were used in operation 104 to determine the generic numerical delay models. Output 110, i.e., the subset of library cells with their specific numerical delay models and optionally the new generic numerical delay models for the library cell types, can then be provided as input to the first phase of numerical synthesis 112.

Examples of systems and techniques that can be used to exclude library cells (which can be used to perform operation 108) are described in U.S. patent application Ser. No. 13/479,807, entitled “Excluding library cells for delay optimization in numerical synthesis,” the contents of which are herein incorporated by reference in their entirety.

The first phase of numerical synthesis 112 sizes cells in circuit design 114 using the subset of library cells with their specific numerical delay models and the new generic numerical delay models. In some embodiments, the first phase of numerical synthesis 112 models a numerical optimization problem based on circuit design 114 and the specific and generic numerical models for the subset of library cells. Output 116 from the first phase of numerical synthesis 112 includes the numerically optimized size for each cell instance in circuit design 114. Specifically, in some embodiments, output 116 includes the numerically optimized (and desired) input capacitance value for each cell instance in the circuit design 114.

Output 116 is then provided as one of the inputs to the second phase of numerical synthesis 118. The other inputs to second phase of numerical synthesis 118 include the library cells in technology library 102 and output 106 which includes the numerical delay models that were generated by operation 104. The second phase of the numerical synthesis 118 then instantiates cells that have the numerically optimized and desired cell sizes in circuit design 114 to produce optimized circuit implementation 120 using the discrete library cells from the technology library. In particular, for each cell instance in circuit design 114, the second phase of numerical synthesis 118 selects a library cell that best matches the numerically optimized size, (if a library cell with exactly the optimized size is not available, a library cell with nearly the optimized size is selected) from the set of library cells that were part of output 106. In some embodiments, second phase of numerical synthesis 118 also receives output 110 as input (shown in FIG. 1 using a dashed line). Specifically, second phase of numerical synthesis 118 uses the fact that certain cells were excluded from output 106 to determine which cells to instantiate. If the user marks certain library cells as “do not use” cells, then the optimization process does not use them.

The reason a subset of library cells is used in the first phase and the set of all library cells is used in the second phase is as follows. The first phase determines numerically optimized cell sizes by solving a numerical optimization problem. In this phase, it is desirable to restrict the search space to library cells that can be accurately modeled using a logical effort and parasitic delay based numerical modeling approach. Therefore, a subset of library cells is used in this phase to restrict the search space. Once the numerically optimized cell sizes have been determined, it is important to enable the circuit synthesis process to select cell instances from a wide range of library cells. Therefore, in the second phase, the restriction on the search space is relaxed and the set of library cells that was in output 106 is used.

In general, gates can be optimized in any arbitrary order. Some embodiments described herein optimize gates in a reverse-levelized order so that, when a current set of gates are being optimized, the loads on those gates can be modeled accurately as fixed loads because those loads would have already been optimized in the current iteration. In a reverse-levelized processing order, each logic gate is associated with a level so that each logic gate's level is greater than the highest level associated with logic gates that are electrically coupled to the logic gate's fan-in. Optimizing gates in reverse-levelized order means that the gates are optimized in decreasing order of their levels, i.e., a gate at level n−1 will only be optimized after all gates at level n have been optimized. Gates at the same level may be optimized in any arbitrary order.

FIG. 2 illustrates a reverse-levelized processing order in accordance with some embodiments described herein. In the figure, the notation “L=n,” where n is a number denotes that the level associated with the logic gate is equal to n. Circuit design 200 includes logic gates 202-214 that have been assigned levels 0 through 3. Logic gates 202-204, whose fan-in is coupled to a primary input of the circuit design, are assigned level 0. All other logic gates are assigned a level that is greater than the highest level associated with the logic gates that are electrically coupled to the logic gate's fan-in. For example, logic gate 206 is assigned level 1, logic gates 208-210 are assigned level 2, and logic gates 212-214 are assigned level 3. In a reverse-levelized optimization order, the optimization process optimizes logic gates in decreasing order of their levels. For example, some embodiments described herein would optimize logic gates 212 and 214 before optimizing logic gates 208 and 210, and so forth.

The above description of the reverse-levelized order was provided for illustration purposes only, and is not intended to limit the embodiments described herein to a particular gate optimization order. It will be apparent to a person having ordinary skill in the art that the gates can be optimized in any desired order.

Modeling the Gate Optimization Problem

FIG. 3A how a gate optimization problem can be modeled in accordance with some embodiments described herein.

The portion of the circuit design shown in FIG. 3A includes gates 302, 304, 306, 308, 310, and 312. The output of driver gate 302 is electrically connected to the inputs of gates 304, 306, and 308. The outputs of gates 304 and 306 are electrically connected to the inputs of gates 310 and 312, respectively. The electrical behavior of the wires that electrically connect outputs and inputs is modeled using a network of resistances and capacitances as shown in FIG. 3A. For example, the electrical behavior of the wire that electrically connects the output of driver gate 302 with the inputs of gates 304, 306, and 308 is modeled using a network of resistances R1-R4 and capacitances CW1-CW4. Likewise the electrical behavior of the wire that electrically connects the output of gate 304 with the input of gate 310 is modeled using resistance R5 and capacitance CW5, and the electrical behavior of the wire that electrically connects the output of gate 306 with the input of gate 312 is modeled using resistance R6 and capacitance CW6.

The input capacitances (i.e., of the inputs shown in FIG. 3A) of gates 302-308 are Cs, C1, C2, and Cƒ, respectively. The labels O₁, O₂, and F₁ correspond to the output pins of gates 304, 306, and 308, respectively. The loads at O₁, O₂, and F₁, are L₁, L₂, and L_(f), respectively. The generic logical effort values for gates 304 and 306 are g₁ and g₂, respectively. The generic parasitic delay values for gates 304 and 306 are p₁ and p₂, respectively. The specific logical effort and parasitic delay values for gate 302 are g_(s). and p_(s), respectively. The specific logical effort and parasitic delay values for gate 308 are g₃ and p₃, respectively.

Delays from pin A to output pins O₁, O₂, and F₁, are modeled as follows. Let the wire delays from pin Z to pins O₁, O₂, and F₁, be dw1, dw2, and dw3. The wire delays dw1, dw2, and dw3 can be modeled using an Elmore delay model. Systems and techniques for efficiently computing Elmore wire delay are described in U.S. patent application Ser. No. 13/485,600, entitled “Incremental Elmore delay calculation,” the contents of which are herein incorporated by reference in their entirety. Specifically, let CWT=CW1+CW2+CW3+CW4. Delay dw1 can be modeled as dw1=τ1+B11·C1+B12·C2, where, τ1=R1·(CWT+Cƒ)+R2·CW2, B11=(R1+R2), and B12=R1.  (4)

Likewise, delays dw2 and dw3 can be modeled as dw2=τ2+B21·C1+B22·C2 and dw3=τ3+B31·C1+B32·C2, respectively, where, τ2=R1·(CWT+Cƒ)+R3·CW3, B21=R1, B22=R1+R3, τ3=R1·(CWT+Cƒ)+R4·(CW4+Cƒ), B31=R1, and B32=R2.  (5)

Let the delays from pin A to output pins O₁, O₂, and F₁, be ƒ₁, ƒ₂, and ƒ₃, respectively. Delay ƒ₁ can be expressed as follows:

$\begin{matrix} \begin{matrix} {f_{1} = {{\frac{g_{s}}{Cs} \cdot \left( {{CWT} + {C\; 1} + {C\; 2} + {Cf}} \right)} + p_{s} + {{dw}\; 1} + {\frac{g_{1}}{C\; 1} \cdot L_{1}} + p_{1}}} \\ {= {{\frac{g_{s}}{Cs} \cdot \left( {{CWT} + {Cf}} \right)} + {\tau 1} + p_{s} + p_{1} + {\left( {\frac{g_{s}}{Cs} + {B\; 11}} \right) \cdot}}} \\ {{C\; 1} + \left( {\frac{g_{s}}{Cs} + {B\; 12}} \right) + {C\; 2} + {\frac{g_{1} \cdot L_{1}}{C\; 1}.}} \end{matrix} & (6) \end{matrix}$

Equation (6) can be rewritten as:

$\begin{matrix} {{f_{1} = {{K\; 1} + {A\;{11 \cdot C}\; 1} + {A\;{12 \cdot C}\; 2} + \frac{\phi 1}{C\; 1}}},} & (7) \end{matrix}$ where,

$\begin{matrix} {{{{K\; 1} = {{\frac{g_{s}}{Cs} \cdot \left( {{CWT} + {Cf}} \right)} + {\tau 1} + p_{s} + p_{1}}},{{A\; 11} = {\frac{g_{s}}{Cs} + {B\; 11}}},{{A\; 12} = {\frac{g_{s}}{Cs} + {B\; 12}}},{and}}{{\phi 1} = {g_{1} \cdot {L_{1}.}}}} & (8) \end{matrix}$

Likewise, ƒ₂ can be expressed as follows:

$\begin{matrix} {{f_{2} = {{K\; 2} + {A\;{21 \cdot C}\; 1} + {A\;{22 \cdot C}\; 2} + \frac{\phi 2}{C\; 2}}},} & (9) \end{matrix}$ where,

$\begin{matrix} {{{{K\; 2} = {{\frac{g_{s}}{Cs} \cdot \left( {{CWT} + {Cf}} \right)} + {\tau 2} + p_{s} + p_{2}}},{{A\; 21} = {\frac{g_{s}}{Cs} + {B\; 21}}},{{A\; 22} = {\frac{g_{s}}{Cs} + {B\; 22}}},{and}}{{{\phi 2} = {g_{2} \cdot L_{2}}},}} & (10) \end{matrix}$ and, ƒ₃ can be expressed as follows: ƒ₃ =K3+A31·C1+A32·C2,  (11) where,

$\begin{matrix} {{{{K\; 3} = {{\frac{g_{s}}{Cs} \cdot \left( {{CWT} + {Cf}} \right)} + {\tau 3} + {\frac{g_{3} \cdot L_{f}}{Cf}p_{s}} + p_{3}}},{{A\; 31} = {\frac{g_{s}}{Cs} + {B\; 31}}},{and}}{{A\; 32} = {\frac{g_{s}}{Cs} + {B\; 32.}}}} & (12) \end{matrix}$

Let D1, D2, and D3 be the desired delays from pin A to pins O₁, O₂, and F₁, respectively. The objective function ƒ_(c) for the gate optimization problem can be defined as follows: ƒ_(c)=max[(ƒ₁ −D1),(ƒ₂ −D2),(ƒ₃ −D3)],  (13) where, “max” is a function that outputs the maximum value from the set of values that are provided as arguments to the function. The term “(ƒ₃−D3)” in Equation (13) can limit the solver from increasing C₁ and C₂. Therefore, in some alternate embodiments, the objective function can be formulated without using the “(ƒ₃−D3)” term. In these embodiments, the capacitance CW4 and the capacitance Cƒ are still used for computing the delay of the driver gate and for computing wire delays from pins Z to pins O₁ and O₂. The following discussion is based on the objective function shown in Equation (13).

If the numerical solver can handle non-differentiable objective functions (such as the “max” function), then the objection function shown in Equation (13) can be provided to the numerical solver to determine optimal gate sizes. Some numerical solvers require that the objective function be differentiable. Specifically, the conjugate-gradient method requires that the objective function be differentiable. Further details of conjugate-gradient based numerical solvers can be found in M. R. Hestenes, and E. Stiefel, “Methods of Conjugate Gradients for Solving Linear Systems,” Journal of Research of the National Bureau of Standards, Vol. 49, No. 6, pp. 409-436, December 1952, and W. C. Naylor, R. Donelly, L. Sha, “Non-linear optimization system and method for wire length and delay optimization for an automatic electronic circuit placer,” U.S. Pat. No. 6,671,859, Dec. 30, 2003, the contents of which are incorporated by reference.

For numerical solvers that require a differentiable objective function, a differentiable function that approximates the behavior of the “max” function in Equation (13) can be used. For example, in some embodiments, the log-sum-exp differentiable function can be used to model the “max” function as follows: ƒ_(c)=α·ln [e ^((ƒ) ¹ ^(−D1)/α) +e ^((ƒ) ² ^(−D2)/α) +e ^((ƒ) ³ ^(−D3)/α)].  (14)

The plot of the objective function shown in Equation (14) has two or more substantially linear regions that have different slopes, and that are connected by curves. The constant α is a scaling factor that affects the smoothness of the curve as the function transitions from one substantially linear region to another.

For large values of α, the connecting curves are wide (i.e., the curves have a large radius of curvature). As α gets smaller, the connecting curves become tighter (i.e., the curves have a small radius of curvature). The function shown in Equation (14) can be made to model the non-differentiable “max” function with an arbitrary level of accuracy by making α arbitrarily small. However, decreasing a can increase the amount of time that the numerical solver requires to converge to a solution. In fact, if α is too small, it can prevent the numerical solver from converging to a solution altogether. Some embodiments described herein use a value of α that is small enough so that the objective function shown in Equation (14) models the “max” function with sufficient accuracy, and is large enough so that the numerical solver is able to converge to a solution within a reasonable amount of time.

Some embodiments use a conjugate-gradient based numerical solver to determine the optimum gate sizes. For example, in FIG. 3A, a conjugate-gradient based numerical solver can be used to determine optimum gate sizes for gates 304 and 306 by determining the optimum input capacitance values C1 and C2, respectively. Specifically, the following expressions can be used to compute the partial derivatives of the objective function ƒ_(c) with respect to capacitances C1 and C2:

$\begin{matrix} {{\frac{\partial f_{c}}{{\partial C}\; 1} = \frac{\mspace{25mu}{{{\mathbb{e}}^{{({f_{1} - {D\; 1}})}/\alpha} \cdot \left( {{A\; 11} - \frac{\phi 1}{\left( {C\; 1} \right)^{2}}} \right)} + {{{\mathbb{e}}^{{({f_{2} - {D\; 2}})}/\alpha} \cdot A}\; 21} + {{{\mathbb{e}}^{{({f_{3} - {D\; 3}})}/\alpha} \cdot A}\; 31}}}{{\mathbb{e}}^{{({f_{1} - {D\; 1}})}/\alpha} + {\mathbb{e}}^{{({f_{2} - {D\; 2}})}/\alpha} + {\mathbb{e}}^{{({f_{3} - {D\; 3}})}/\alpha}}},{and}} & (15) \\ {\frac{\partial f_{c}}{{\partial C}\; 2} = {\frac{\mspace{25mu}{{{{\mathbb{e}}^{{({f_{1} - {D\; 1}})}/\alpha} \cdot A}\; 12} + {{\mathbb{e}}^{{({f_{2} - {D\; 2}})}/\alpha} \cdot \left( {{A\; 22} - \frac{\phi 2}{\left( {C\; 2} \right)^{2}}} \right)} + {{{\mathbb{e}}^{{({f_{3} - {D\; 3}})}/\alpha} \cdot A}\; 32}}}{{\mathbb{e}}^{{({f_{1} - {D\; 1}})}/\alpha} + {\mathbb{e}}^{{({f_{2} - {D\; 2}})}/\alpha} + {\mathbb{e}}^{{({f_{3} - {D\; 3}})}/\alpha}}.}} & (16) \end{matrix}$

A conjugate-gradient based numerical solver iteratively optimizes the objective function. In each iteration of the numerical solver, the current values of capacitances C1 and C2 are used to determine the partial derivatives (e.g., by using Equations (15)-(16)), and the partial derivatives are used to determine a new set of values for capacitances C1 and C2. The numerical solver outputs a solution when a termination condition is reached. For example, the numerical solver may output a solution when the relative and/or absolute change in the solution vector (e.g., the capacitances C1 and C2) is less than a pre-determined threshold. It will be apparent to a person having ordinary skill in the art that other termination conditions (e.g., a time limit, an iteration limit, etc.) can also be used either separately or in concert. Note that the partial derivatives can be computed incrementally, e.g., by computing the values of certain terms in Equations (15)-(16) only once and then reusing the values of those terms in subsequent calls to the conjugate-gradient solver.

The numerical solver can converge faster if a seed that is close to the optimal solution is provided to the numerical solver. Some embodiments use a closed analytical expression to estimate optimal sizes of the input capacitances of the gates that are being optimized (e.g., capacitances C1 and C2 in FIG. 3A), and provide the estimated input capacitance values as a seed (i.e., a starting point) to the numerical solver. Systems and techniques for estimating gate sizes are described in U.S. patent application Ser. No. 13/537,880, entitled “Estimating optimal gate sizes by using numerical delay models,” the contents of which are herein incorporated by reference in their entirety.

Although some embodiments described herein used the Elmore wire load model, the systems and techniques described in this disclosure can generally use any wire load model. Specifically, the systems and techniques described herein can use a best case wire load model (which is a model that ignores the wire resistance and models the wire load using a single lumped capacitance), a worst case wire load model (which is a model that uses a lumped resistance and a lumped capacitance to model the wire load), and a balanced-tree wire load model (which is a model that distributes the total capacitance and resistance of the wire equally across each branch, and uses a lumped capacitance and lumped resistance to model the wire load of each branch). FIGS. 3B-3D illustrate a best case wire load model, a worst case wire load model, and a balanced-tree wire load model, respectively, in accordance with some embodiments described herein. In FIGS. 3B-3D, RW is the wire resistance and the CW is the wire capacitance.

FIG. 4 illustrates a process for optimizing a circuit design in accordance with some embodiments described herein. An embodiment (e.g., a computer system) can perform the process, which begins by receiving a circuit design (operation 402). The circuit design can generally be at any abstraction level, e.g., the circuit design can be a logical design or a physical design.

Next, the embodiment begins processing gates in the circuit design in a given processing order (operation 404). For example, in some embodiments, the embodiment can process the gates in a reverse-levelized processing order as explained in reference to FIG. 2. The embodiment then selects a portion of the circuit design for optimization (operation 406). For example, the embodiment may select the portion of the circuit design shown in FIG. 3A for optimization.

Next, the embodiment can model a gate optimization problem for the selected portion of the circuit design (operation 408). Modeling the gate optimization problem can involve collecting (e.g., by retrieving values stored in memory locations) circuit information for the selected portion of the circuit design, and constructing expressions for determining values of entities that are required during optimization.

In some embodiments, the selected portion of the circuit design can include a driver gate that drives one or more inputs of each gate in a set of gates. In these embodiments, the collected circuit information can include generic logical effort values of each gate in the set of gates, an input capacitance value and a specific logical effort value of the driver gate, and wire resistance and capacitance values of a net that electrically connects an output of the driver gate with one or more inputs of each gate in the set of gates. The circuit information can further include a generic logical effort value for a virtual driver gate, and a fixed load value as seen by the virtual driver gate.

Embodiments that use a conjugate-gradient based numerical solver may construct expressions for computing the conjugate gradient (e.g., Equations (15)-(16)). Note that if the selected portion of the circuit design includes n gates that are being optimized simultaneously, then the embodiments compute n partial derivatives with respect to the n input capacitances corresponding to the n gates that are being optimized. Some or all parts of the “collecting information” operation and/or the “constructing expressions” operation can be performed once and then reused later if required.

The embodiment can then determine a seed for the numerical solver by estimating optimal gate sizes using a set of closed-form expressions (operation 410). Operation 410 is optional, i.e., some embodiments may directly proceed to operation 412 without performing operation 410. A seed is a starting point in the solution space. For example, if the numerical solver is simultaneously optimizing n gate sizes, then the seed could be a set of n input capacitances corresponding to the n gate sizes that are being optimized. In general, the closer the seed is to the optimal solution, the faster the numerical solver is expected to converge to the optimal solution. Some embodiments determine the seed by using the systems and techniques for estimating gate sizes that are described in U.S. patent application Ser. No. 13/537,880, entitled “Estimating optimal gate sizes by using numerical delay models.”

Next, the embodiment solves the gate optimization problem using a numerical solver (operation 412). For example, some embodiments that use a conjugate-gradient based numerical solver can perform operation 412 as follows. The embodiments can determine the conjugate-gradient based on the collected circuit information. The conjugate gradient is defined by the set of partial derivatives of the objective function (e.g., Equation (14)) with respect to the variables (e.g., input capacitances C1 and C2 of gates 304 and 306, respectively) that are being optimized. Next, the embodiment can provide the value of the conjugate gradients as well as the value of the objective function to the numerical solver, and the numerical solver can return with a new set of values for the variables (e.g., input capacitances C1 and C2) that are being optimized. The embodiment can then compute the new conjugate gradients, and provide the new conjugate gradients and the value of the objective function to the numerical solver. This iterative search for the optimal solution can terminate when the difference in the objective function between successive solutions returned by the numerical solver is less than a given threshold.

Once the selected portion of the circuit design has been optimized, the embodiment can check if all gates have been processed (operation 414). If more gates need to be processed (“No” branch), then the embodiment can select the next portion of the circuit design for optimization (operation 416) and return to operation 408.

On the other hand, if all gates have been processed (“Yes” branch), then the embodiment can check if more iterations are required (operation 418). If no more iterations are required (“No” branch, e.g., because the gate sizes have converged to a substantially optimal value or some other termination condition has been met), then the embodiment can output the optimized circuit design (operation 420) that includes gates with the optimized gate sizes. On the other hand, if more iterations are required (“Yes” branch), the embodiment can return to operation 404 and begin processing gates in the circuit design in a given processing order (with the new gate sizes that were determined in the previous iteration).

FIG. 5 illustrates a computer system in accordance with an embodiment of the present invention.

A computer system can generally be any system that can perform computations. Specifically, a computer system can be a microprocessor, an application specific integrated circuit, a distributed computing system, a cloud computing system, or any other computing system now known or later developed. Computer system 502 comprises processor 504, memory 506, and storage 508.

Computer system 502 can be coupled with display 514, keyboard 510, and pointing device 512. Storage 508 can generally be any device that can store data. Specifically, a storage device can be a magnetic, an optical, or a magneto-optical storage device, or it can be based on flash memory and/or battery-backed up memory. Storage 508 can store application 516, operating system 518, and data 520.

Application 516 can include instructions that when executed by computer system 502 cause computer system 502 to perform one or more processes that are implicitly or explicitly described in this disclosure. Data 520 can include any data that is inputted into or outputted by application 516.

The above description is presented to enable any person skilled in the art to make and use the embodiments. Various modifications to the disclosed embodiments will be readily apparent to those skilled in the art, and the general principles defined herein are applicable to other embodiments and applications without departing from the spirit and scope of the present disclosure. Thus, the present invention is not limited to the embodiments shown, but is to be accorded the widest scope consistent with the principles and features disclosed herein.

The data structures and code described in this disclosure can be partially or fully stored on a non-transitory computer-readable storage medium and/or a hardware module and/or hardware apparatus. A non-transitory computer-readable storage medium includes all computer-readable storage mediums with the sole exception of a propagating electromagnetic wave or signal. Specifically, a non-transitory computer-readable storage medium includes, but is not limited to, volatile memory, non-volatile memory, magnetic and optical storage devices such as disk drives, magnetic tape, CDs (compact discs), DVDs (digital versatile discs or digital video discs), or other media, now known or later developed, that are capable of storing code and/or data. Hardware modules or apparatuses described in this disclosure include, but are not limited to, application-specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), dedicated or shared processors, and/or other hardware modules or apparatuses now known or later developed.

The methods and processes described in this disclosure can be partially or fully embodied as code and/or data stored in a non-transitory computer-readable storage medium or device, so that when a computer system reads and executes the code and/or data, the computer system performs the associated methods and processes. The methods and processes can also be partially or fully embodied in hardware modules or apparatuses. Note that the methods and processes can be embodied using a combination of code, data, and hardware modules or apparatuses.

The foregoing descriptions of embodiments of the present invention have been presented only for purposes of illustration and description. They are not intended to be exhaustive or to limit the present invention to the forms disclosed. Accordingly, many modifications and variations will be apparent to practitioners skilled in the art. Additionally, the above disclosure is not intended to limit the present invention. The scope of the present invention is defined by the appended claims 

What is claimed is:
 1. A method for optimizing a circuit design, the method comprising: optimizing, in a computer system, a portion of the circuit design that includes a driver gate that drives one or more inputs of each gate in a set of gates, wherein said optimizing the portion of the circuit design comprises: modeling a gate optimization problem for the portion of the circuit design based on circuit information that includes generic logical effort values of each gate in the set of gates, an input capacitance value and a specific logical effort value of the driver gate, and wire resistance and capacitance values of a net that electrically connects an output of the driver gate with one or more inputs of each gate in the set of gates, wherein said modeling includes constructing a differentiable objective function that minimizes a maximum delay from an input of the driver gate to an output of each gate in the set of gates; executing a conjugate-gradient based numerical solver to solve the gate optimization problem; incrementally computing a gradient of the differentiable objective function; and providing the computed gradient to the conjugate-gradient based numerical solver.
 2. The method of claim 1, further comprising: estimating gate sizes for the set of gates based on a set of closed-form expressions, and providing the estimated gate sizes as a seed to the conjugate-gradient based numerical solver to solve the gate optimization problem.
 3. The method of claim 1, wherein the wire resistance and capacitance values of the net are based on an Elmore wire load model.
 4. The method of claim 1, wherein the wire resistance and capacitance values of the net are based on a best case wire load model.
 5. The method of claim 1, wherein the wire resistance and capacitance values of the net are based on a worst case wire load model.
 6. The method of claim 1, wherein the wire resistance and capacitance values of the net are based on a balanced-tree wire load model.
 7. A non-transitory computer-readable storage medium storing instructions of an electronic design automation (EDA) software tool that, when executed by a computer system, cause the computer system to perform a method for optimizing a circuit design, the method comprising: optimizing a portion of the circuit design that includes a driver gate that drives one or more inputs of each gate in a set of gates, wherein said optimizing the portion of the circuit design comprises: modeling a gate optimization problem for the portion of the circuit design based on circuit information that includes generic logical effort values of each gate in the set of gates, an input capacitance value and a specific logical effort value of the driver gate, and wire resistance and capacitance values of a net that electrically connects an output of the driver gate with one or more inputs of each gate in the set of gates, wherein said modeling includes constructing a differentiable objective function that minimizes a maximum delay from an input of the driver gate to an output of each gate in the set of gates; executing a conjugate-gradient based numerical solver to solve the gate optimization problem; incrementally computing a gradient of the differentiable objective function; and providing the computed gradient to the conjugate-gradient based numerical solver.
 8. The non-transitory computer-readable storage medium of claim 7, wherein the method further comprises: estimating gate sizes for the set of gates based on a set of closed-form expressions, and providing the estimated gate sizes as a seed to the conjugate-gradient based numerical solver to solve the gate optimization problem.
 9. The non-transitory computer-readable storage medium of claim 7, wherein the wire resistance and capacitance values of the net are based on an Elmore wire load model.
 10. The non-transitory computer-readable storage medium of claim 7, wherein the wire resistance and capacitance values of the net are based on a best case wire load model.
 11. The non-transitory computer-readable storage medium of claim 7, wherein the wire resistance and capacitance values of the net are based on a worst case wire load model.
 12. The non-transitory computer-readable storage medium of claim 7, wherein the wire resistance and capacitance values of the net are based on a balanced-tree wire load model.
 13. A computer system, comprising: a processor; and a storage medium storing instructions of an electronic design automation (EDA) software tool that, when executed by the processor, cause the computer system to perform a method for optimizing a circuit design, the method comprising optimizing a portion of the circuit design that includes a driver gate that drives one or more inputs of each gate in a set of gates, wherein said optimizing the portion of the circuit design comprises: modeling the gate optimization problem for the portion of the circuit design based on circuit information that includes generic logical effort values of each gate in the set of gates, an input capacitance value and a specific logical effort value of the driver gate, and wire resistance and capacitance values of a net that electrically connects an output of the driver gate with one or more inputs of each gate in the set of gates, wherein said modeling includes constructing a differentiable objective function that minimizes a maximum delay from an input of the driver gate to an output of each gate in the set of gates; executing a conjugate-gradient based numerical solver to solve the gate optimization problem; incrementally computing a gradient of the differentiable objective function; and providing the compound gradient to the conjugate-gradient based numerical solver.
 14. The computer system of claim 13, wherein the method further comprises: estimating gate sizes for the set of gates based on a set of closed-form expressions, and providing the estimated gate sizes as a seed to the conjugate-gradient based numerical solver to solve the gate optimization problem.
 15. The computer system of claim 13, wherein the wire resistance and capacitance values of the net are based on a wire load model that is one of: an Elmore wire load model; a best case wire load model; a worst case wire load model; and a balanced-tree wire load model. 